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Abstract 

We present a calculation of the neutrino-nucleon scattering cross section which takes into 
account the nuclear correlations in the relativistic random phase approximation. Our approach 
is based on a quantum hadrodynamics model with exchange of a, uj, it, p and 5 mesons. In view 
of applications to neutrino transport in the final stages of supernova explosion and protoneutron 
star cooling, we study the evolution of the neutrino mean free path as a function of density, 
proton-neutron asymmetry and temperature. Special attention was paid to the issues of renor- 
malization of the Dirac sea, residual interactions in the tensor channel and meson mixing. It 
is found that RPA corrections, with respect to the mean field approximation, amount to only 
10% to 15% at high density. 

1 Introduction 

Neutrino transport is a key ingredient to understand the mechanism of supernova explosion and 
the subsequent cooling of the proto neutron star formed in the collapse. In the delayed explosion 
mechanism suggested by Wilson, the problem of shock stall can be overcome if sufficient energy 
is deposited by the neutrinos to revive the shock. It is seen in recent numerical simulations ^, 
H|, [I], ||, H of core collapse that, besides other issues such as convection, the final outcome of the 
explosion depends sensitively on the rate of neutrino energy deposition. 

The following stage of protoneutron star cooling also depends on the rate of neutrino trans- 
port. It determines the shape of the neutrino signal emitted during the first few seconds of the 
protoneutron star formation. It has also been suggested that the deleptonization could trigger an 
instability as the star cools down, and lead to the delayed formation of a black hole, in order to 
explain the time interval observed between the 8th and 9th neutrino detected in the explosion of the 
SN1987a supernova. It is thus important to refine the models, since we now begin to benefit of the 
observation of young neutron star surface temperatures, and of the upgrading of neutrino detectors 
such as SuperKamiokande, MACRO, SNO or AMANDA (see e.g. R) to collect observational data 
from an eventual supernova explosion. 

An enhanced neutrino emission can occur at high densities if the neutrino-nucleon cross section 
is modified through the nuclear correlations, as was recognized in the early works of Iwamoto and 
Pethick H and Sawyer Q. Interest in this topic has recently been revived by the feasability of full 
Boltzmann simulations of neutrino transport in collapsing supernovae j|, |3|, |||, ||] . 

While early calculations only were taking into account the nucleon mass reduction in the 
medium, important progress has been realized since then in order to include Hartree-Fock fl(| , 



random phase approximation (RPA) 0, ||, ||, 0, |T|, ||, |7| or ladder ||, ||, Jx± |2l], ||, ||, || 
correlations. It is generally found that the neutrino opacities are suppressed by medium effects, 
if however no collective mode is excited [16]. As a matter of fact, the latter authors find a size- 
able enhancement due to the excitation of a spin zero sound mode, announcing a transition to a 
ferromagnetic state in a nonrelativistic treatment. Due to the degree of uncertainty existing on 
the equation of state of neutron star matter at high density, the existence of such a transition 
remains an open question, although it seems unfavoured in a relativistic treatment. Even if no 
such instability exists, it is important to check the order of magnitude of the suppression factor to 
be applied to the scattering rate. 

In order to asses the importance of some points raised in previous calculations, we present in 
this paper a new calculation of the neutrino-nucleon rate in dense, hot, asymmetric matter. After 
a short presentation of the relativistic meson exchange model and relevant formulae in section 
two, we first perform a few calculations in a simplified model where the background matter is 
approximated by pure neutron matter. In particular, we examine the effect of the Dirac sea and 
of various renormalization procedures commonly used in the litterature in order to subtract the 
divergences in the vacuum fluctuation term. As a matter of fact, calculations reported in [25] 
adscribed a sizeable correction due to this effect. In contrast, we found here only a minor effect on 
the size and position of the zero sound mode in the longitudinal component, which moreover gets 
washed out by Landau damping and integration on the exchanged momentum and energy. These 
results are presented in section three. 

We also studied the way of introducing the short-ranged residual interaction in the tensor 
channel, which has to be taken into account in order to avoid the appearance of unobserved pion 
condensation at low density. In the non relativistic formalism, this is commonly done by adding a 
Landau-Migdal contact interaction. In the relativistic case, this has been done by various methods, 
one beeing Horowitz's Ansatz in the pion propagator [12, 26]. Two previous works [13, 17] using 
this Ansatz found their results to be very sensitive to the choice of the numerical value of the 
Landau-Migdal constant g'. We compare results obtained with this Ansatz to those obtained from 
a more transparent and straightforward way of introducing this term through a contact Lagrangian. 
We conclude that there is no dependence on the value of g' , once the contact term is properly taken 
into account, and point out a possible explanation for the discrepancy with the authors previously 
quoted. These results are presented in section four. 

In a second part of this paper, corresponding to section four, we study the full model where 
the background matter is asymmetric. The proton fraction may be given as a fixed value, or be 
determined by (3 equilibrium in the two extreme situations of neutrino free matter Y u = 0, or 
trapped neutrinos with a lepton fraction typical of supernovae Yl ~ 0.4. Here we point out the 
fact that the phenomenon of meson mixing has not been taken fully into account in the previous 
treatments. Indeed, in addition to the usual a-io mixing, the neutron-proton asymmetry also opens 
the possibility of a-p and UJ-p mixing in the case of neutral current relevant for neutrino-nucleon 
scattering. Finally, we consider here the contribution of the 5 meson. This meson is usually left 
out, since the a, to and p mesons are sufficient to describe the properties of nuclear matter at the 
mean field level. In asymmetric matter however, it can have a non negligible contribution and it 
mixes with the a, u and p mesons. 
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2 Neutrino-nucleon scattering in the relativistic random phase 
approximation 

2.1 Free scattering rate 

Two related processes are needed in order to calculate the mean free path and spectrum of neutrinos 
in matter. They are the neutrino-nucleon scattering through the neutral current 



u{K) + n(P) -> v(K') + n(P') 
and the absorption (or creation by the inverse process) through the charged current 

v(K)+n{P)^e-{K')+p(P') 
In either case, the neutrino-nucleon cross section may be written as 

da 



G F E v i 



dE u dn 64vr 3 E v 



(1) 



(2) 



(3) 



We will concentrate in this work on the neutral current process. 

Note that we do not include in this definition the Pauli blocking factor for the outgoing lepton. It 
will be taken into account in the final stage of the calculation when we compute the mean free path 
of the neutrino in matter. In this formula, the tensor represents the lepton current, coupled 
with the weak vertex T' We = 7 a (l — 75). For neutral current scattering with massless neutrinos, we 
have 



L a/3 



Tr 



We 



2K a K (B + (K.q)g a(B - (K a qP + q a K< i ) T ie a/3 ^ u K^q v 



(4) 



K and K' are the four momenta of the ingoing and outgoing leptons. For neutral current scattering 
we have K 2 = K' 2 = 0. The scattering angle is defined by K.K' = E v E' v cos9. E v is the ingoing 
and E v i = E v — uj is the outgoing lepton energy. The energy loss to is the zero component of 
momentum exchange = — K 1 ^. 

The structure function can be related to the imaginary part of the retarded polarization 



S^(q) = Jd 4 xe iqx < J"{x)r(0) > 
-2 



1 - e" 

and the differential cross section takes the form 

da 



G 2 F E v i 



1 



dE v dtt 32vr 3 E v 1 - e" 



(5) 
(6) 

(7) 



The factor (1 — e 



with z = — Afi) arises from detailed balance. Afi is the difference 



between the chemical potential of the outgoing and ingoing nucleons. For the neutral current 
process, A/i = 0. For the charged current process, Afi = (x = fi n — fi p . At the mean field level, H^f 
is the (retarded) polarization 



ImUf 



tanh 



n 



a/3 
11 



:Jd 4 pTi [T a G u (p)T /3 G 11 ( P + q) 



(8) 
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with T a being the weak vertex to the hadronic current T a = j a (Cv — C/iTs)- For the neutral 
current, C% = -1/2, C\ = -g A /2 for the neutron and C v = 1/2 - 2 sin 2 6 W , C A = 9a/2 for the 
proton, with g A = 1.23 the axial coupling constant and sin 2 0\y = 0.232 is the Weinberg angle. For 
the charged current we have Cy = cos# c , Ca = gACOs6 c , with cos# c = 0.95 the Cabbibo angle. 
G n (p) is the nucleon propagator. At the mean field level, the nucleon behaves as a quasiparticle 
with effective mass M, momentum P and chemical potential /i, and the propagator can written as 



the sum of a vacuum and a density dependent term (see e.g. [27, 28]) 

G l \p) = ( 7 .P + M) | p2 _ l M2 - - + 2m8(P 2 - M 2 ) [9(p Q )n(p Q ) + fl(-po)rc(po)]} (9) 

with n ^°) = e PivX) + l ' "(PP) = e -/3( P0 -,) + 1 ( 10 ) 
2.2 RPA corrections 

In the hot and dense medium through which the neutrino is moving in the protoneutron star, the 
nucleon is not free, but correlated with the other nucleons through the strong interaction. Several 
effects contribute, among which are: 

— The fact that the nucleons can be described in the medium as quasiparticles with effective masses 
and chemical potentials. This represents the mean field approximation. This correction has already 
been taken into account in the preceding section in the definition of the nucleon propagator. 

— One can go a step further and include the exchange term at the level of Hartree-Fock correlations. 



This is the approach of e.g. Fabbri and Matera [1C]. We will not consider this type of corrections 
here. 

— Another consequence of the correlations is the broadening of the nucleon width, corresponding 



to the inclusion of ladder corrections fl8|, 19, 2C, E2J. This is the equivalent for scattering to the 



Landau- Migdal-Pomeranchuk effect studied in the case of neutrino Bremsstrahlung by [g^, 24 , |2^] . 
— Finally, one has to take into account the RPA correlations. This is the subject of the present 
work. 

The strong interaction may be modelized by relativistic a, u>, tt and p meson exchange. We will 



also include the 5 meson in a later stage of this calculation (see section §4.1) in order to investigate 
some specific effects related to the fact that the matter of the proto neutron star has a sizeable 
proton-neutron asymmetry. We will work in a relativistic formalism with an interaction Lagrangian 
of the quantum hadrodynamics type: 

Ant = V> [~9crO + Qujl^fi ~ ^-in^d^.T - g s S.T + gp'fpn-r + ~^ alJV ^Pn-^j i> 

--bm N (g a a) 3 - ^c{g a a) 4 (11) 

The coupling of the pion is taken in the pseudovector form, since this is known to reproduce better 
the phenomenology of nucleon-pion scattering. It will be seen anyway, at the end of the calculation, 
that in fact the pion does not contribute directly to the neutral current process. The non-linear 
a couplings <t 3 , <r 4 are introduced in order to obtain a better description of the equation of state 
with a value of the incompressibility modulus and effective mass at saturation density in agreement 
with the experimental data. 

RPA correlations can be introduced by substituting the mean field polarization in Eq. (0) by 
the solution of the Dyson equation 

LL WW — LL WW LL WS u SSa(3 LV SW 

a,b=cr,uj,p,7r... 
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where the index W stands for a vertex with a weak coupling and S for a vertex with a strong 
coupling. The tilde II indicates a resummed polarization. D$s is the propagator of the mesons 
a = a, u>, p°, 7T° for the neutral current process (we will see later that the pion does not contribute) 
or a = p^, ir^ for the charged current process. The first term of the Dyson equation corresponds 
to the mean field approximation taken in the preceding section. 
The Dyson equation can be rewritten as 

LL WW — LL WW LL WS u SSaf3 LV SW \ L6 ) 

a,b=<r,uj,p,7r... 

which provides a solution in terms of the meson propagator D$s dressed in the RPA approximation. 
2.3 Polarizations 



Polarization insertions were calculated in the relativistic random phase approximation. They are 
given by the loop integral 



n 



AB 



d A p 



Tr 



T A G(p)T B G(p + 



(14) 



G(p) is the nucleon propagator calculated in the Hartree approximation, as given in Eq. (|). The 
F A are the vertices of the interaction, which can be either of the weak type T\y = CyJ^ + CaIsI* 1 
or of the strong type T$ with S = {<r, u, p, ir, 5}. Formulae are given for the polarizations in e.g. 
1| 111 || 1§ |§. 

The polarizations can also be obtained by performing a linear response analysis on a Hartree 
ground state. This was the method used in the following references, to which the reader is referred 
for a full account of the calculation techniques and meson dispersion relations. The method was 
described in [40|, with results given for the ir meson in [p4[ |, for the a, and u> mesons in [35, 36] 
and for the p in [37, 38]. The notation we are using here can be found in [39|, which also presents 
a discussion of some results concerning the introduction of the contact term. Results concerning 
asymmetric matter are presented in |40[|. A new feature, which has been overlooked in previous 
determinations of the RPA corrections [|13|, |i~7"| , [OJ, is the mixing between mesons of different isospin 
when the distribution functions of the protons and neutrons differ. As a consequence, the Dyson 
equation ([D^) acquires a more complex structure. The formulae needed for the calculation of the 
mixed polarizations and meson propagators were given in [ficfl . 

The polarizations calculated in the linear response analysis coincide with those available in 
the litterature from the Green function formalism at zero temperature and in symmetric nuclear 
matter. At finite temperature, the real parts coincide, but some precisions are necessary concerning 



the imaginary parts. In the real time formalism [28], one defines the Green function, and self 
energies as 2x2 matrices with the indices labelling the two branches of the time contour. The 
various components are related to each other. In this paper, we chose to work with the retarded 
polarization. 

Polarizations are generally obtained under the form of a "density dependent" and a "vac- 
uum fluctuation" term. This latter part entails divergences which have to be subtracted by some 
renormalization procedure. Regularization of divergences ocurring in vacuum polarization may 
be performed through the inclusion of a counterterm lagrangian with couplings adjusted so as to 
cancel the infinities. The subtraction is defined up to some residual finite constant, which is then 
determined by choosing a renormalization point where the polarization and its derivatives with 
respect to the momentum and effective mass is required to vanish. 
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The criteria guiding this choice will have to be reexamined in the light of the recent advances 
concerning effective field theories. As a matter of fact, the various possible renormalization schemes 
unfortunately lead to discrepancies observed among the results available in the litterature. A 
particularly distressing example was the case of the p meson discussed in J38[ . One cannot ignore 
the problem, since simply eliminating the vacuum by a normal reordering produces pathologies in 
the dispersion relation. The problem roots in the fact that we are working in an effective theory, 
which is non renormalizable in the usual QED sense. A renormalization can still be performed at 
a given order. Nevertheless, the parameters of the theory, which are adjusted so as to reproduce 
available experimental data, should enforce the symmetries and scalings of the integrated degrees 
of freedom of the underlying more fundamental theory. There is some hope that one could solve 
the problem using "naturalness" arguments | |41j] . Fortunately, as we will see in section §||, our final 
result has only a negligibly weak dependence in the choice of the renormalization condition. 

The polarization which enters the definition of the differential neutrino-nucleon scattering cross 
section may be decomposed onto orthogonal projectors, formed with the vectors and tensor available 
in the problem, i.e. the metric g^ v = diag (1, — 1, — 1, — 1), the hydrodynamic velocity and the 
transferred momentum (f 3 . 

U^ w = fl T + U L + U Q + i U E (15) 



with 

rpflU 



rftf 
rj 2 

LIU _ 



if 



rfrf 



e^Vpqx 



(16) 



The imaginary unit i appearing before the antisymmetric tensor E^ combines with the i present 
in the corresponding term of the lepton current, so that the product involves the imaginary part 
of the polarization II^;. The following properties might be useful 



A^A^ : 



: T^k c 

A"" , 

: E^ 



= A^Qa = A^E^ = E>* a Q» = 
E» a E a u = q 2 7] 2 T^ 



(17) 



When calculating the polarization of the mesons with formula (14), a zero sound branch ap- 
pears in the dispersion relation of the pion. Even though it is weaker than in the nonrelativistic 
formulation, it remains a spurious effect, since no such effect is observed experimentally. When 
analyzing the structure of the pion potential, it is seen that this is related to the short range be- 
havior. The introduction of the short range effects actually should come out of a full many-body 
calculation. Nevertheless, it is a common practice to implement it by the definition of a residual 
interaction. In the non relativistic formalism, one adds the Landau-Migdal contact interaction 
Vc = g' {f-w/m- w ) 2 5{r)tJi.02. In the relativistic case, several Ansatze have been suggested. For 



example, Horowitz el al. [12, 26 1 modify the free pion propagator as follows 

1 



mi 



mi 



and obtain the pion and rho propagators dressed at RPA level by inverting G 1 
explicitely 



G, 



-i 



rr 



QflU 

QP- 



-1 



G 





n {Q)p,v 



P 



m v 



i -i 



(18) 



i.e. 



(19) 
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In these expressions, the polarizations H% p and Ii^ u are calculated with a pseudowector "pion" 
vertex (fn/mw)^!' 1 - The true pion is obtained by projecting over q^. 

Another more conventional possibility is to add a contact term to the Lagrangian 

£ 3 <?' (Jf) 2 (^757^) (^757^) (20) 

When the dispersion relations of the mesons are calculated using the method of linear reponse 
analysis on the Hartree ground state, this leads to self consistent equations for the pion and rho 
polarizations |39| whose solution is given by (pi, 22). 

Both alternatives modify the polarization of the pion and of the transverse mode of the p meson 
as follows 

- -A- (21 > 

q 2 - g>Il n 

n pT — n pT - vq9 ._ pv (22) 



1 + g'IL^ 

with U nT = --U^T lxv and rfq 2 ^ = --R^iE. 



However, the first procedure presents the drawback of introducing additional terms in the 
dispersion relation related to the unphysical components of the auxiliary "pion". These terms 
would be responsible for the appearance of spurious branches in the spacelike region. A more 
detailed account of this discussion can be found in |39! 



3 Preliminary study in pure neutron matter 

The RPA correction to the mean field value of the polarization entering the definition of the 
neutrino-nucleon differential scattering cross section is 



ah^pa 



n 



(a)H jj(w)/LtQ jjip)^ TT-M^t 



ws 



ws 



n 



ws 



n 



ws 



/ 







o \ 























V o 











1L sw 

lL SW 
\ LL SW 



AU L RPA A^ + AU T RPA T^ + AU Q RPA Q^ + iAU E RPA E^ 



(23) 



Introducing the development of the polarizations and propagators on the projection tensors ([16]) 
we obtain (where we have dropped the "WS" index indicating that the polarization loops connect 
a weak vertex with a strong one) 



An£ PA 

^RPA 



An 
An 



RPA 



RPA 



U a G a Uar] 2 + 2 n^G^n^ L rj 2 -U^lG^ lU^ L -U pL G p L U P L (24) 

—U w tGu T^iuj T + ~RweG uj Tn EV 2 q 2 - UpTGprllpT + U pE G pT U pE rj 2 q 2 (25) 

— IH^tGuT^-ujE — tUpTGpTTlpE (26) 

-n^n, (27) 



An^ pA must now be contracted with the tensor of the lepton current L pu . Using the property 



(28) 



we can see that the pion does not contribute to the neutrino-nucleon scattering. 
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3.1 Response functions 

The contraction of the lepton current with the polarization can be expressed by means of three 
structure functions Ri(q), R 2 {q) and Rs(q) defined as follows |jT^l 

S^(q) = J d 4 xe iq - x < J^(x)r(0)> (29) 
= R 1 (q)uV + R 2 (q)(uV - gT) + # 3 (?)<f<f + R^u" + <fu v ) + iR 5 (q)e^ r, a q^ 
Using the properties 

L 00 = 8E V {E U - uj)(cos6 + 1) , = 16K.q = 16E u E' u (cos9 - 1) (30) 

L^K„ V = -j^-L 00 = -^8E„E' v (co S 6 + 1) (31) 

8uj 2 + 8E U E' U (3 - cos 0)1 (32) 
E^L^ = 8q 2 (E v + E' u ) (33) 



ZTTr = 8q 2 + = g 



we arrive at 

1m (L^n« 



1 - e 



4E U E U , [#i(l + cos 6) + i? 2 (3 - cos 9) - 2{E V + E U ,)R 5 (1 - cos 9)] (34) 



The structure functions are related to the previous polarizations IP" = II^p + An^ PA 

= n L + n T t» u + n Q + » n B by 



R\ = 



1 


- e- z 




2 


1 


- e" 2 




2 


1 


- e" 2 



<? 2 w 2 „ 



(35) 



i? 2 = -Im[n T ] (36) 

1 — e 2 

i? 5 = — Im[n E ] (37) 
1 — e 2 

In the non relativistic limit, R± and -R 2 reduce to the density and spin density correlation functions 
respectively. The axial-vector structure function R§ appears only in a relativistic treatment. 

The "vector" structure function i?i involves the sigma meson and the longitudinal part of the 
omega meson, as well as the longitudinal part of the rho. The contribution of the latter meson 
is negligible in this channel. As is known from previous studies, a zero sound mode appears in 
the mixed ct-uj dispersion relation. It manifests itself as a pole in the RPA meson propagator. At 
finite temperature, it is quenched by Landau damping. As an example, the longitudinal response 
function was represented on Fig. 1 at twice the saturation density and at a temperature T = 30 
MeV for a fixed exchanged 3-momentum k = 100 MeV. The longitudinal response function is 
suppressed in the RPA approximation (full line) as compared to the mean field approximation 
(dashed line). The zero-sound appears as a peak at the high frequency edge of the distribution 
strength. Note, however, that this mode will lie mostly outside of the integration range when we 
calculate the mean free path. Indeed, the fact that the exchanged momentum must be spacelike 
restricts the frequency to uj 2 < k 2 ; moreover, the condition on the scattering angle — 1 < cos# < 1 
imposes that k 2 < (2E V — uj) 2 . The integration range is therefore restricted to k £ [0,oo[ and 
uj £ [— k, min(fc, 2E U — k)] or, equivalently, uj G [—oo,E u ] and k G [|u;|,2ii^ — uj]. 

The RPA correction to the transverse polarization is mostly determined by the p meson, whereas 
the a and uj mesons contribution to this quantity remains small. An example of the shape of the 
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Fig. 1 - Structure function 
Rl, displaying overall RPA re- 
duction and zero sound en- 
hancement 
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transverse structure function is displayed on Fig. 2 in the mean field approximation (full line) and 
with RPA corrections included (dashed line). The results of Fig. 2 were obtained for a density 
twice that of saturation, a temperature 30 MeV and at a fixed value of the 3-momentum 100 MeV. 

We find that the dependence in the Landau-Migdal parameter g' is very small, when it is 
introduced through a contact Lagrangian (CL), in contrast with the results reported in [13, 17]. 
The dashed line representing our RPA result is unchanged at the level of precision of the figure 
whether we take g' = 0.6 or g' = 0. A tentative explanation for this discrepancy could be the fact 
that the later references [O, 17] are using Horowitz Ansatz (HA). Indeed, when the w-p mixings 



and transversal projection of the Yl% u tensor introduced through the inversion of Eq. ( |19|) are kept 
in the equations, we would obtain the dotted line of Fig. 2 represented for g' = 0.6 These terms 
actually would dominate over the other contributions from the p meson, and would be responsible 
for a strong suppression of the transverse response function. They are proportional to g'. As was 
argued in the preceding section and in [39] however, these contributions are spurious and should 
be subtracted. This could be the reason for the strong g' dependence observed by the authors of 
p"3| , p"7| ], and of the large reduction factors obtained in the RPA approximation in these references. 
We believe that this point needs to be investigated further. 

As discussed in the previous section, the polarizations entering the definition of the dressed 
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Fig. 3 - Structure function 
R5, displaying a shift towards 
lower frequencies 



co [MeV] 



meson propagators in the RPA approximation contain a contribution from the vacuum fluctua- 
tions which needs to be renormalized. At the present time, there exist some discrepancies in the 
litterature concerning the choice of a renormalization procedure. As discussed e.g. in |39[ , 
this is a source of uncertainty in the determination of the dispersion relations of the mesons in the 
framework of hadronic models, and in particular in the prediction of the behavior of the effective 

We refer the reader to 



p meson mass in the medium. We refer the reader to |38|, [4C| |4^, 43 1 for details. Before we 
proceed, it is necessary to know in what measure does this issue affect the results presented in 
this work. Accordingly, we calculated the response functions in three cases. A first case consists 
of simply dropping the vacuum term (which is equivalent to perform a normal ordering). In the 
second case, ("scheme A" in this work) we chose a renormalization scheme as used by Kurasawa 
and Suzuki [E^j which minimizes the coupling to counterterms on the q 2 = shell. In the third 
case, we used a renormalization procedure which preserved the formal structure of the expressions 
of the vacuum polarizations as a function of the effective nucleon mass ("scheme 3" of [38], called 
"scheme B" in this work). 
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Fig. 4 - Effect of renormal- 
ization scheme on longitudinal 
structure function 



The choice of the renormalization procedure has some small influence on the longitudinal R\ 
response function. As seen on Fig. 4, its main effect is to modify the strength and position of the 



10 



zero-sound mode. It was however argued above that the contribution of the zero sound mode to the 
neutrino-nucleon scattering is only marginal, since its characteristic frequency lies mostly outside 
of the values kinematically allowed for neutrino-nucleon scattering. Moreover, as will be seen in the 
next paragraph, the longitudinal reponse is not the dominant contribution to the neutrino-nucleon 
scattering. The structure functions R2 and R5 are not appreciably affected by the choice of the 
renormalization scheme. The difference between one or the other schemes are of the order of the 
fraction of a percent and therefore not distinguishable on the scale of the figures. As a consequence, 
we are able to conclude that the vacuum fluctuations will not affect the neutrino-nucleon scattering 
opacity in any appreciable way. 



3.2 Differential cross section 

The differential cross section can now be calculated from Eq. (0) by inserting the previous results 
for the longitudinal, transverse and axial vector polarizations. It is instructive to compare the 
relative magnitude of their contributions, which are displayed in Fig. 5. The result is rescaled by 
dividing it by the textbook estimate o"o = G 2 F (C v + 3C\) E^/ir Q]. The transverse contribution is 
dominant and will provide for about 60 % of the total result. We find that it is very little modified 
by RPA correlations. The corrections arise from the subdominant longitudinal and axial-vector 
polarizations ITx and He- The longitudinal contribution is affected by a reduction factor of about 
two. We notice at the right and left ends of the figure two small peaks which arise from the zero- 
sound mode excitation. The contribution of the axial-vector response function is only marginal at 
the density (ns = 2 n sat ) chosen in this figure, although it is larger at higher density. It is slightly 
enhanced by RPA correlations and its strength is shifted towards lower values of the energy loss 
of the neutrino uj = E u — E' u , thus acting to counterbalance the reduction from the longitudinal 
contribution. 



T=30 MeV 

n B= 2n Sc 

k=25 MeV 



MF-L 
RPA-L 

MF-T 
RPA-T 



-0.8 




0.8 



Fig. 5 - Relative contributions 
of longitudinal, transverse and 
axial vector polarizations to the 
differential scattering cross sec- 
tion 



As a conclusion to this preliminary study, we gather here the main results of the section. We are 
able to report that the choice of renormalization scheme does not affect the results more than on 
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the level of less than 1 %. Further good news is that our results are also insensitive to the choice of 
the residual contact (Landau-Migdal) interaction. The transverse response function is left almost 
unchanged by RPA correlations, while the longidudinal and axial vector responses are reduced by 
about a factor of two. Since the transverse structure function yields the dominant contribution to 
the neutrino-nucleon scattering, we do not expect spectacular modifications from RPA correlations 
as compared to the mean field result. 



4 Full model - Asymmetric nuclear matter 

Up to now, we have made the simplifying assumption of pure neutron matter. We now pass to 
describe a more realistic model, where the proton fraction is defined by the physical conditions in 
the star. The proton fraction is defined by 

Y p = (38) 

Pp + Pn 

where p p and p n are the proton and neutron densities. Numerical calculations show that the 
chemical equilibrium is reached very rapidly. The proton fraction will therefore be determined 
by (5 equilibrium. As a rule of thumb, in cooled neutron stars, the neutrinos can leave the star 
unhindered, and we have Y p ~ 0.1. On the other hand, in supernovae the neutrinos are still trapped 
dynamically inside the matter on the diffusion time scale and contribute to displace the equilibrium 
to higher ratios of the proton fraction Y p ~ 0.3. 

The asymmetry enters at several levels in the calculation of the neutrino nucleon scattering 
cross section. It first enters in the determination of the thermodynamics: we have kp p ^ kp n , and 
also M p 7^ M n in a model with 5 meson exchange. It also enters in the polarizations, e.g we must 
sum the proton and neutron contributions in polarizations such IT- = U^ 1 ^ +11^ . Moreover, new 
mixing channels occur which are not present in symmetric matter. They arise from polarizations 
involving vertices with mesons of different isospin, which are given by the difference between the 
proton and neutron contributions, such as e.g. Ii ap = Il^f^ — n^p"^ in longitudinal modes or 
II^p = nlfp ' — Tl^p^ in transverse modes. 



4.1 Thermodynamics 

The basic properties of nuclear matter can be reproduced at the mean field level with a, uj and p 
meson exchange, when non linear sigma meson couplings l/36m<r 3 + l/4ccr 4 are taken into account 
in order to obtain a better value of the compressibility (see e.g. j45fl ). 

The 5 meson also appears in meson exchange models such as the Bonn potential |46|, fiTp . Since 
the 5 meson carries isospin, it can give important contributions in strongly asymmetric matter. 
Moreover, when density dependent mean field models are adjusted in order to reproduce the most 
recent Dirac-Brueckner-Hartree-Fock calculations, it is seen that a 5 meson is needed in order to 
reproduce the results at finite asymmetry H8L E9[. The behavior of the equation of state when a 



5 meson is present was studied by Kubis and Kutschera 5Cj. The 5 meson is at the origin of a 



difference between the neutron and the proton effective masses. 

M n = m- g a o + g$6 , M p = m - g a o - g 5 5 (39) 

On the other hand, it brings a negative contribution to the asymmetry energy, so that the coupling 
to the p meson has to be readjusted to a higher value to compensate for that effect. 

In this work, the strength of the 5 coupling g$ will be fixed to the Bonn potential value |4(| 
/47r = 1.1075. The couplings g a , g w , g p , b, c were adjusted [40] so as to reproduce the experimental 
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saturation point with 

po = 0.17 fm~ 3 , B/A = -16MeV, K = 250 MeV, m*/m = 0.8, a A = 30 MeV (40) 
and take the values 

g a = 8.00, <? w = 7.667, 6 = 9.637 10~ 3 , c = 7.847 10~ 3 , <? p = 4.59 (41) 

The ratio k p = f p /g p was fixed to the Bonn value n p = 6.1. 

The proton fraction is determined by the f3 equilibrium condition jl = \i n — fi p = fi e — [i v . 
If neutrinos are trapped inside matter, the chemical potential of the neutrino has a finite value 
\i v = (67r 2 pYu) 1 /^ . For a typical value of the lepton fraction Yl = Y v + Y e = 0.4, [i v is of the 
order of 200-250 MeV and the proton fraction Y p of the order of ~ 0.3 - 0.36. Yl is determined by 
neutrino transport (e.g., diffusion equation), to which the neutrino-nucleon cross section serves as 
input. If the matter is transparent to neutrinos, fi v = and the proton fraction is of the order of 
~ 0.1. The following figures (Fig. 6) show the behavior of thermodynamical parameters relevant 
for the calculation of the neutrino opacities. 
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The neutron star matter is assumed to be in (3 equilibrium. The two cases of neutrino free 
matter (Y u = 0) or matter with trapped neutrinos (Yj, = 0.4) are displayed. All figures were drawn 
assuming a vanishing temperature. The upper left panel shows the evolution of the proton fraction 
as a function of the baryon density. The upper right panel shows the behavior of the proton and 
neutron effective masses. The difference M p — M n is larger in neutrino free matter since it allows for 
a smaller proton fraction and therefore a larger mean 5 field. On the lower panel are represented the 
difference between the chemical potentials of the neutron and the proton £i in the cases (Y u = 0) 
and (Yl = 0.4) as a full line. In matter with a trapped lepton fraction Yl = 0.4, the neutrino 
chemical potential is non vanishing; it is compared to p on Fig. 6-(c). The full line was obtained 
for a coupling to the 5 field gs = 5. If this coupling is set to zero, the chemical potentials are given 
by the dashed lines on Fig. 6-(c). 



4.2 Meson mixing 

The dispersion relations of the a-u)-5°-p° system of neutral mesons, which is relevant for the calcu- 



lation of the neutrino-nucleon scattering through the neutral current process, were studied in 
As discussed in this reference, the difference existing between the neutron and proton distribution 
functions makes mixing possible in all chanels in asymmetric matter. The pion does not mix with 
the other mesons, and, as pointed out in section HQ its contributions to the u-N scattering process 
vanishes, so we will not consider it further. The RPA correction is now obtained from 
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Explicit expressions for the dressed meson propagator matrix were given in [40|. The calculation 
proceeds as in section §||. The longitudinal projection of the RPA correction to the polarization 
now receives contributions from the 5 meson and from all possible combinations of the mixings, 
that is, not only from a-u and 5-p already present in symmetric matter, but also from 5-a, a-p, uj-p 
and oj-5 mixings. The transversal T and axial vector E components receive additional contributions 
from uj-p mixing. Finally, the projection parallel to q^q u vanishes through contraction with the 
(massless) lepton current, it would contain the contribution of the pion. 

The expressions of Ira (U^ 1 ' L^ u ) so obtained are introduced in the formula for the differential 
cross section (|7|). We calculated separately the longitudinal, transverse and axial- vector contribu- 
tions to the differential cross section and represent them in Fig. 7 as a function of the ratio of 
the neutrino loss energy to the transferred three momentum k = 50 MeV. The thermodynamical 
conditions were chosen to be those of neutrino free matter (Y u = 0) with a density n = 4n sa t and 
finite temperature T=20 MeV. The transferred 4-momentum is subject to kinematical constraints 
which restrict the range over which it will be integrated in order to obtain the total cross section 
and the mean free path (see next section). The limit of the integration range is also represented 
on Fig. 7. 

It is seen from these figures that the contribution of the longitudinal polarization is appreciably 
affected by meson mixing. The longitudinal modes involve the full a — uj — 5 — p mixing. The 
mixing is stronger at high asymmetry, high density and high momentum transfer. There exists as 
before a zero sound mode, but it is undiscernible on the figures, because it is partially washed out 
by Landau damping and lies outside of the integration range. 

The transverse and axial vector contributions involve only p — uj mixing in the transverse modes. 
The modification to the transverse contribution is small in the kinematically allowed region. On 
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Fig. 7 - Contribution of longitudinal 
(upper left panel), transverse (upper right 
panel) and axial-vector (lower panel) po- 
larizations to the differential scattering 
cross section, showing the correction aris- 
ing from fully taking into account meson 
mixing. 



the other hand the axial vector response function is somewhat enhanced by p-tu mixing. 

The conclusions that were already presented from the preliminary study are confirmed in this 
more elaborate model. We obtain that the dominant contribution to the neutrino-nucleon scattering 
cross section, that is the transverse one, is hardly modified by the RPA corrections as compared 
to the mean field approximation. The RPA corrections affect the longitudinal and axial vector 
components, the longitudinal part beeing suppressed and the axial part enhanced with respect to 
the mean field approximation, and partially cancelling each other. The effect of meson mixing is of 
the order of 20 % on both the longitudinal and axial vector parts, but further act counterbalance 
each other. The net result is that we do not expect dramatic corrections to the global result from 
RPA correlations. 



4.3 Mean free path 

Finally, we will calculate the contribution of neutrino-nucleon scattering to the mean free path, 
since it gives an immediate feeling of the phenomenon of opacity when it is compared with the 
length scales typical of the star, e.g. the star radius or a typical convection length scale. We 
must naturally stress the fact that this is only an estimate. As a matter a fact, several other 
processes, like neutrino emision and absorption through charged current, scattering on electrons, 
Bremsstrahlung ... contribute to the definition of the actual mean free path. 

The mean free path is defined here as the inverse of the total cross section per unit volume 
obtained by integrating the differential cross section calculated in the preceding sections. 

1 - ° 2p 1 r ada P"" ±s (1 ~ fiK)) Zm ( L^U R ) US) 
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with 



mm(2E u — k, k) 



(1 — f{E' v ) = (1 + exp[{E' u — fi v )/T])) is a Pauli blocking factor for the outgoing lepton. The 
chemical potential \i v is determined by (3 equilibrium and by neutrino transport equations {e.g., 
diffusion equation), which determine the lepton fraction Yj, at a given instant of the protoneutron 
star cooling. 

In the following figure (Fig. 8), we represented the ratio of total cross sections as calculated in 
the RPA and mean field approximations. The chosen thermodynamical conditions are representa- 
tive of the earlier stage of the cooling, when the neutrinos are still trapped inside the matter (the 
lefpton fraction was chosen to be Yl=0A, for typical values of the temperature (T = 20 MeV) and 
density (n = 4n sat ). The neutrino energy was fixed to E u = 30 MeV. 
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At high density, the total neutrino-neutron scattering cross section is found to be reduced by 
RPA correlations by a factor (only) ~ 10%. When full meson mixing is taken into account, the 
reduction is even smaller. This value should be compared to suppression factors of the order of 
two quoted in recent publications 13, 17]. A tentative explanation for this discrepancy was given 
in section §j||. Further investigation on this issue is needed. 

The mean free path obtained from the integration of fl43|) is shown on Fig. 9 in matter with 
trapped neutrinos as a function of density for various values of the temperature. At high density, 
the mean free path is somewhat lenghtened by RPA correlations as compared to the mean field 
value. At low density and moderate temperature, on the other hand, RPA correlations would yield 
an enhancement of the cross section and a reduction of the mean free path. A similar behavior 



was reported in [17]. It should be kept in mind, however, that the validity of the model becomes 
questionable in this range. 



5 Summary and perspectives 

Let us gather here the main results obtained in this work. Special attention was paid to the issues 
of renormalization of the Dirac sea, residual interactions in the tensor channel and meson mixing. It 
has been shown that the vacuum fluctuations have only a negligible effect on the scattering rate. We 
examined two different prescriptions to introduce the contact term in the tensor channel, and argue 
that, when properly taken into account, the numerical value of the Landau- Migdal parameter does 
not affect the structure functions. As a consequence, the transverse contribution to the scattering 
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rate is basically unchanged, while the contributions coming from the longitudinal and axial-vector 
parts combine to produce a 10% to 15 % reduction of the scattering rate with respect to the mean 
field approximation. 

A natural extension of this work concerns the absorption and emission cross section by the 
charged current reaction in the RPA approximation. This is presently under study. 

As the density increases, a more complete description should include the hyperons A, 
Taking into account these baryons does not present any particular difficulty, since the structure of 
the equations is the same, only with different values of the coupling constants. This has been done 
in [33, 51 1 . As these authors argue moreover, the high fraction of trapped neutrinos and relatively 
hot temperatures both conspire to strongly suppress the formation of exotic, strangeness carrying 
particles, so that these are not expected to have a strong influence on the first few seconds of the 
life of the protoneutron star. Nevertheless, the exoticas will be determinant on the later stage of 
the protoneutron star cooling, at t ~ 30 — 60 s, which are marginally or will soon become observable 
with the next generation of neutrino detectors now on construction, such as UNO. 

A more important problem is to arrive to an accurate description of the short range correlations. 
These have been studied separately from RPA corrections by several autors (see eg. [19, 2(], 24]). 
The calculated correction factors to the neutrino opacities are large, especially in the low density 
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regime. This is in fact related to the well-known problem of nuclear matter descriptions, which 
should simultaneously take into account short and long range corrections, "ladders" and "loops", 
in a consistent way. While methods do exist to do so fl52| , at least in an approximate way, they are 
quite unwieldy, especially in the relativistic formulation. Some results in this direction have been 
reported recently in . 

In matter at subnuclear density, and also at high density, the homogeneous phase may not 
be the most stable state of matter. It will therefore be important to study the role of ordered 
configurations, as coherent scattering is expected to dominate in this regime. At low density, and 
sufficiently low temperatures, the spinodal instasbility triggers condensation of droplets of a denser 
phase in a more dilute gas principally composed of neutrons. The droplets take a spatially ordered 
configuration due to Coulomb forces. This will eventually form the crust of the cooled neutron 
star. At high density, there are various possibilities of forming ordered structures. For example, if 
a transition to quark gluon plasma takes place, an important fraction of the matter of star can be 
in a mixed phase, and will acquire an ordered structure similar to that present in the crust fl45|| . 
Another possibility is the formation of a pion condensate in the alternating spin layer configuration 



53]. In any case, such structures are generally unfavored by temperature and high lepton (and 
proton) fractions, so that the same remark as done before applies: this type of correction are 
expected to set in at a later phase of the protoneutron star cooling. It could be very interesting 
once we are able to detect the tail of the neutrino emission from a supernova event; if a sudden 
change in the neutrino would occur, it could be interpreted as a signal that a phase transition has 
taken place. 
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